Non-Markovian homodyne-mediated feedback on a two-level atom: a quantum 

trajectory treatment 
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Quantum feedback can stabilize a two-level atom against decoiierence (spontaneous emission), 
putting it into an arbitrary (specified) pure state. This requires perfect homodyne detection of the 
atomic emission, and instantaneous feedback. Inefficient detection was considered previously by two 
of us. Here we allow for a non-zero delay time r in the feedback circuit. Because a two-level atom is a 
nonlinear optical system, an analytical solution is not possible. However, quantum trajectories allow 
a simple numerical simulation of the resulting non-Markovian process. We find the effect of the time 
delay to be qualitatively similar to that of inefficient detection. The solution of the non-Markovian 
quantum trajectory will not remain fixed, so that the time-averaged state will be mixed, not pure. 
In the case where one tries to stabilize the atom in the excited state, an approximate analytical 
solution to the quantum trajectory is possible. The result, that the purity {P — 2Tr[p^] — 1) of the 
average state is given by P = 1 — 47r (where 7 is the spontaneous emission rate) is found to agree 
very well with the numerical results. 

42.50.Lc, 42.50.Ct, 03.65.Bz 



I. INTRODUCTION 

The theory of open quantum systems is fundamental to 
understanding dissipation on a microscopic and macro- 
scopic level in diverse fields (measurement theory, quan- 
tum optics, quantum chaos, solid state, quantum compu- 
tation) [|l]-^, wherever quantum irreversibility matters. 
Traditionally, open quantum systems are described by 
the reduced density operator, which is obtained from the 
total density operator by tracing over the environmental 
degrees of freedom. The dynamics of the reduced density 
operator is commonly described by a master equation. 

Besides the density matrix formalism, there has been 
an increasing interest over the last years in the quantum 
trajectory approach which uses a stochastic Schrodinger 
equation for the state vector |-0) |6|-||]. This approach 
serves as a numerical tool for solving the master equation, 
as the reduced density operator is recovered as the ensem- 
ble average of these stochastically evolving pure states. 
It also has a deeper significance, as, in some cases, the 
quantum trajectories can be interpreted as the evolution 
of the system state conditioned on continuous monitoring 
of its environment . This means that quantum tra- 
jectories are ideal for treating quantum feedback, where 
the measurement results from this monitoring are used 
to control the dynamics of the system 

The history of feedback-control in open quantum sys- 
tems goes back to 1980's with the work of Yamamoto 
and co-workers [|l^,|l^, and Shapiro and co-workers [ p^ . 
Their objective was to explain the observation |l^] of sub- 
shot-noise fluctuations in an in-loop photocurrent. They 
did this using quantum Langevin equations (stochastic 
Heisenberg equations for the system operators) and also 
semiclassical techniques. The latter approach was made 



fully quantum- mechanical by Plimak For systems 

with linear dynamics, all of these approaches, and the 
quantum trajectory approach of Refs. |ll|,|l^, are equally 
easy to use to find analytical solutions. The advantage 
of the Wiseman and Milburn (quantum trajectory) ap- 
proach to quantum control via feedback is for systems 
with nonlinear dynamics, as we will discuss later. 

Before proceeding further, it is necessary to distinguish 
what we mean by 'quantum control' via feedback from an 
alternative usage of the term 'quantum control' which has 
become current in laser chemistry |]l8| . The work in this 
field has sought sufficient control over laser pulses so that 
very specific quantum dynamics of electronic and vibra- 
tional states may be achieved, primarily with the objec- 
tive of making and breaking bonds in molecular systems. 
Typically such experiments are conducted on a large en- 
semble of identical prepared target systems. Two kinds 
of control are distinguished; (i) learning control, and, (ii) 
feedback control. In learning control the objective is to 
design the required laser pulse as follows. After apply- 
ing a pulse to a sample of many identical constituents 
(an ensemble), one measures some property of the sys- 
tem, adjusts the laser pulse according to the measure- 
ment results, then applies the new laser pulse to a new 
sample and so on. The cycle of adjust-excite-measure is 
repeated until the desired results are obtained. In feed- 
back control, the laser excitation always acts on the same 
sample of many identical systems without re-preparation, 
through the same repeated sequence. Of course the ad- 
justment one makes to the next pulse will be quite dif- 
ferent in the two cases. 

In contrast to this controlling of the average state of 
an ensemble, we are concerned with the control of the 
state of a single quantum system subject to continuous 
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measurement. While this is similar to (ii), the fact that 
only a single quantum system is involved, not an ensem- 
ble, makes it fundamentally different. The quantum effi- 
ciency of the measurement has to be relatively large for 
the control to be at all effective. This means that one 
has to take into account the quantum back-action of the 
measurement on the system, which is not necessary in 
the laser chemistry 'feedback control' referred to above. 
Only very recently has it become possible to isolate, and 
continuously interrogate a single quantum system at the 
back-action limit. This has been achieved in a number 
of fields including cavity QED ion trapping |^ and 
mesoscopic electronics 

Returning to the different approaches to control of 
a single quantum system, another advantage of the 
Wiseman-Milburn quantum feedback theory is that it is 
very easy to consider the limit of Markovian (i.e. instan- 
taneous) feedback. In this limit it is possible to derive a 
master equation which describes the unconditioned sys- 
tem dynamics including the effect of feedback QJl^ . In 
previous work, two of us applied Markovian quan- 
tum feedback theory to a two-level atom, where the light 
emitted by the atom was detected and used to control the 
dynamics of the atom. Using the feedback master equa- 
tion, analytical results are easily derived. We showed 
that, with one exception, the feedback can generate an 
arbitrary stable pure state of the atom. This is an ex- 
ample of using feedback to control decoherence (due to 
spontaneous emission). 

Decoherence in open quantum systems refers to the 
tendency for pure system states to become entangled 
with many different states of the environment [p3| . If 
the states of the environment are then averaged out, the 
system state tends to become mixed. The radiative decay 
of a two level atom with a non-zero dipole is an example. 
In this case the atomic state becomes correlated with the 
electromagnetic field states. 

Our previous work was motivated by the 1998 work 
of Hofmann, Mahler, and Hess (HMH) Q, in which it 
was shown that by making part of the coherent driving 
of a two-level atom proportional to the homodyne pho- 
tocurrent, it was possible to stabilize the state at any 
point on the bottom half of the Bloch sphere. Our work 
p2| , using the Wiseman-Milburn feedback theory, rean- 
alyzed their proposal and generalized their results in two 
ways. First, we showed that any point on the upper or 
lower half (but not the equator) of the Bloch sphere may 
be stabilized. Second, we considered non-unit-efficiency 
detection, and quantified the effectiveness of the feedback 
by calculating the maximal purity obtainable in any par- 
ticular direction in Bloch space. 

It is now of interest to extend our theory of feedback- 
mediated decoherence control to a non-Markovian pro- 
cess. This would arise naturally in an experiment, since 
any feedback apparatus would necessarily have some de- 
lay associated with it. It is also of interest theoretically 
because non-Markovian feedback in a nonlinear system 
has not been considered in detail before. The two-level 



atom, being the simplest possible system with nonlinear 
dynamics, is an ideal testing ground for descriptions of 
non-Markovian processes. Of course the master equa- 
tions derived in Refs. cannot be applied to non- 
Markovian feedback in which the delay is significant on 
the time scale of the system. 

Most of previous approaches to non-Markovian feed- 
back have been limited to linear systems such as an op- 
tical cavity mode with feedback based on homodyne- 
detection. Wiseman and Milburn used a quantum 
trajectory treatment, as did Dohcrty and Jacobs [ p6[ . 
Giovannetti, Tombesi and Vitali p?! used the quantum 
Langevin approach of input-output theory. As noted 
above, linear systems allow analytical solutions, and all 
approaches to quantum feedback work equally as well. 
However, for a nonlinear system like a two-level atom 
an analytical solution for non-Markovian feedback is not 
possible in general. Hence a numerical solution is nec- 
essary. The only practical way of doing such a simula- 
tion is to use the stochastic quantum trajectories for the 
conditioned system state, which underly the Wiseman- 
Milburn feedback theory. Such a numerical quantum 
trajectory simulation was done recently for a non-linear 
system consisting of a particle in a double-well potential 
with feedback determined by LQG (linearized quadratic 
Gaussian) control theory ||2q] . 

Quantum trajectories describing non-Markovian feed- 
back are of some intrinsic interest as they are a special 
example of non-Markovian quantum trajectories. These 
have been the subject of much interest lately. An ex- 
ample of this is the stochastic wave function method 
developed by Jack, CoUett and Walls where the 

non-Markovian stochastic equation of motion for the 
state vector involves a multiple time integration over 
the system's history conditioned on the measurement 
record over a finite time interval. Furthermore, it has 
been shown recently by Diosi and Strunz and co-workers 
pO| , pl| that it is in principle possible to construct a 
stochastic Schrodinger equation which describes the non- 
Markovian time evolution of any open quantum system. 

Using non-Markovian quantum trajectory simulations 
we are able to determine the effect of a time delay on the 
feedback stabilization scheme introduced in Refs. |p^ , p2| . 
We find that the effect of a time delay is qualitatively sim- 
ilar to that of inefficient detectors. That is, the average 
state can no longer be made pure. States in the upper 
half of the Bloch sphere are more affected than those in 
the lower half, and states near the equator are affected 
most of all. In one special case, where the feedback aims 
to stabilize the atom in the excited state, it is possible 
to linearize the quantum trajectory if the time delay is 
short, and so obtain an approximate analytical solution. 
The result compares very well with the numerical solu- 
tion. It is not obvious that this linearization could be 
performed using any of the other approaches to quantum 
feedback. 

The paper is organized as following. In Sec. II we intro- 
duce the measurement scheme (homodyne detection) for 
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our system (the two- level atom). In Sec. Ill, we briefly 
review our previous work of Markovian feedback stabi- 
lization. The effectiveness of non-Markovian feedback is 
studied in Sec. IV, both numerically and analytically. In 
Sec. V we conclude with a discussion of our results. 



II. THE MEASUREMENT SCHEME 

A. The driven, damped atom 

Consider an atom, with two relevant levels {|5),|e)} 
and lowering operator a = |(7)(e|. Let the decay rate 
be 7, and let it be driven by a resonant classical driv- 
ing field with Rabi frequency 2a. This is as shown in 
Fig. |l|, where for the moment we are omitting feedback 
by setting A = 0. This system is well- approximated by 
the master equation 



(2.1) 



where the Lindblad superoperator is defined as usual 
V[A]B = ABA'' ~ {A^A,B}/2. In this master equation 
we have chosen to define the a-x — cr+a^ and ay — ia—icr^ 
quadratures of the atomic dipole relative to the driving 
field. The effect of driving is to rotate the atom in Bloch 
space around the y-axis. The state of the atom in Bloch 
space is described by the three-vector {x,y,z). It is re- 
lated to the state matrix p by 
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(/ + Xax + yCTy + ZCFz) . 



(2.2) 




FIG. 1. Diagram of the experimental apparatus. The laser 
beam is split to produce both the local oscillator j3 and the 
field Qo which is modulated using the current I{t). The mod- 
ulated beam, with amplitude proportional to a-|-A7(t), drives 
an atom at the center of the parabolic mirror. The fiuores- 
cence thus collected is subject to homodyne detection using 
the local oscillator, and gives rise to the homodyne photocur- 
rent I{t). 

It is easy to s how that the stationary solution of the 
master equation (2.1) is 



-4q;7 



Vss = 0, 



(2.3) 
(2.4) 
(2.5) 



For 7 fixed, this is a family of solutions parameterized by 
the driving strength a G (— (X),cx)). All members of the 
family are in the x-z plane on the Bloch sphere. Thus 
for this purpose we can reparametrize the relevant states 
using r and 9 by 



r sm ( 



r cost 



where 9 e [— tt, tt]. Since 



Tr[p2 



1 



(l + .T^ +l/ + Z^) 



(2.6) 



(2.7) 



is a measure of the purity of the Bloch sphere, r — 
yx^+z^, the distance from the center of the sphere, is 
also a measure of purity. Pure states correspond to r = 1 
and maximally mixed states to r = 0. 



N 




-0.5 



FIG. 2. Locus of the solutions to the Bloch equations. The 
ellipse in the lower half plane is the locus for the equations 
with driving only. The full circle (minus the points on the 
equator) is the locus for the equations with optimal driving 
and feedback, as defined in Sec. III. 
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The locus of solutions in this plane (an ellipse) is shown 
in Fig. ^. Since Zss < 0, all solutions are in the lower 
half of the Bloch sphere. That is, we are restricted to 
\9\ > tt/2. Also, it is evident that the smaller 16*1 is (that 
is, the more excited the atom is), the smaller r is (that 
is, the less pure the atom is). At \9\ = tt the stationary 
state is pure, but this is not surprising as it is simply the 
ground state of the atom with no driving. As \6\ 7r/2 
we have r — > 0. This can only be approached asymptoti- 
cally as |a| — > oo. In summary, the stationary states we 
can reach by driving the atom are limited, and generally 
far from pure. 



B. Homodyne detection 

Now consider subjecting the atom to homodyne de- 
tection. As shown in Fig. |^, we assume that all of the 
fluorescence of the atom is collected and turned into a 
beam (represented in Fig. || by placing the atom at the 
focus of a mirror). Ignoring the vacuum fluctuations in 
the field, the annihilation operator for this beam is ^/7cr, 
normalized so that the mean intensity j (^a^a'^ is equal 
to the number of photons per unit time in the beam. 
This beam then enters one port of a 50:50 beam split- 
ter, while a strong local oscillator (3 enters the other. To 
ensure that this local oscillator has a fixed phase relation- 
ship with the driving laser used in the measurement, it 
would be natural to utilize the same coherent light field 
source in the driving process and as the local oscillator 
in the homodyne detection. This is as shown in Fig. |l|. 

Again ignoring vacuum fluctuations, the two field op- 
erators exiting the beam splitter, bi and 62, are 



h = [Vi^ - (-i)'-/3] /V2 



(2.8) 



When these two fields are detected, the two photocur- 
rents produced have means 

h = {Wf - {-if (V7/3^^ + Vl'^P*) + T^tV) /2. (2.9) 

The middle two terms represent the interference between 
the system and t he l ocal oscillator. 

Equation Eq. (2.9) gives only the mean photocurrent. 



In an individual run of the experiment for a system, what 
is recorded is not the mean photocurrent, but the in- 
stantaneous photocurrent. This photocurrent will vary 
stochastically from one run to the next, because of the ir- 
reducible randomness in the quantum measurement pro- 
cess. This randomness is not just noise, however. It is 
correlated with the evolution of the system and thus tells 
the experimenter something about the state of the sys- 
tem. In fact, if the detection efficiency is perfect, the 
system is collapsed into a pure state, rather than the 
mixed state which is the solution of the master equation. 
The stochastic evolution of the state of the system con- 
ditioned on the measurement record is called a "quan- 
tum trajectory" M. Of course, the master equation is 



still obeyed on average, so the set of possible quantum 
trajectories is called an unraveling of the master equa- 
tion It is the conditioning of the system state on 
the photocurrent record that allows feedback of the pho- 
tocurrent to control the system state. The application 
of an appropriate feedback loop to this continuous mea- 
surement process (to be considered in Sec. Ill) realizes 
an effective "reservoir engineering" to control the system 
at the quantum level. 

The ideal limit of homodyne detection is when the lo- 
cal oscillator amplitude goes to infinity, which in prac- 
tical terms means |/3p ^ 7. In this limit, the rate of 
the photodetections goes to infinity and thus it should 
be possible to change the point process of photocounts 
into a continuous photocurrent with white noise. Also, 
the only relevant quantity is the difference between the 
two photocurrents. Suitably normalized, this is ||6|,p3| 



lit) 



hit) -hit) 



(2.10) 



A number of aspects of Eq. ( ^.10 ) need to be explained. 
First, $ = arg/3, the phase of the local oscillator (de- 
fined relative to the driving field). Second, the subscript 
c means conditioned and refers to the fact that if one 
is making a homodyne measurement then this yields in- 
formation about the system. Hence, any system aver- 
ages will be conditioned on the previous photocurrent 
record. Third, the final term ^(i) represents Gaussian 
white noise, so that 



^it)dt = dWit), 

an infinitesimal Wiener increment defined by 

dWit)^ = dt, 
B[dWit)] = 0. 



(2.11) 



(2.12) 
(2.13) 



Since the stationary solution of the master equation 
confines the state to the x-z plane, it makes sense to 
follow HMH by setting <f> = 0. In that case. 



/(^)-^/7(OcW + eW• 



(2.14) 



That is, the deterministic part of the homodyne pho- 
tocurrent is proportional to Xc = {<^x)c- This should be 
useful for controlling the dynamics of the state in the 
x-z plane by feedback, as we will consider in Sec. III. 
Of course, all that really matters here is the relationship 
between the driving phase and the local oscillator phase, 
not the absolute phase of either. 

The conditioning process referred to above can be 
made explicit by calculating how the system state 
changes in response to the measured photocurrent. As- 
suming that the state at some point in time is pure (which 
will tend to happen because of the conditioning anyway), 
its future evolution can be described by the stochastic 
Schrodinger equation SSE. i,p| 
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d\A{t)) - A,{t)\ip,{t))dt + B,{t)\^P,{t))dW{t). (2.15) 

This is an Ito stochastic equation [Q with a drift term 
and a diffusion term. The operator for the drift term is 



while that for the diffusion is 



(2.16) 
(2.17) 



Both of these operators are conditioned in that they de- 
pend on the system average 



(2.18) 



As stated above , on average the system still obeys the 
master equation (2.1). This is easiest to see from the 
stochastic master equation (SME), which allows for im- 
pure initial conditions. The SME can be derived from 
the SSE by constructing 

d(|7Ae)(V'c|) =(d|^c))(^c| + |^c)(rf(V'c|) 

+ (dm) {d{^,\) , (2.19) 



using the Ito rule ( 2.12| ), and then identifying \^pc){''Pc 
with pc- The result is 



dpc = dtjVlajpc - idta[ay,pc] 



dW{t)y^H[a]pc, 

(2.20) 



where n[A]B = AB + BA^ - Tt[AB + BA^B. Ahhough 
this has been derived assuming pure initial conditions, 
it is valid for any initial conditions [^ . This is also an 
Ito equation, which means the evolution for the ensemble 
average state matrix 



p{t) = E[pe(t)] 



(2.21) 



is found simply by avera ging over the photocurrent noise 
term by using Eq. ( 2.13). T his procedure yields the orig- 
inal master equation (2.1) again. The term "quantum 
trajectory" can be applied to any stochastic conditioned 
evolution of the system, be it described by a SSE or SME. 



III. MARKOVIAN FEEDBACK STABILIZATION 

We now include feedback onto the amplitude of the 
driving on the atom, proportional to the homodync pho- 
tocurrent, as done by HMH. This is as shown in Fig. |^, 
where the driving field passes through an electro-optic 
amplitude modulator controlled by the photocurrent, 
yielding a field proportional to a -I- XI {t). This means 
that the feedback can be described by the Hamiltonian 



In this section we are assuming instantaneous feedback, 
while feedback with time delay will be considered in the 
next section. 

Since the homodync photocurrent ( |2.1C| ) is defined in 
terms of system averages and the noise dW{t), the SSE 
including feedback can still be written as an equation of 
the form ( 2.15| ). The effect of the feedback Hamiltonian 
can be shown (l^j2^ to change the drift and diffusion 
operators to 

Mt) = I [-'tV + (a,), (<)a - {a^)l (0/4 

+ AV7 {^x), {t)ay - 2aV] - AV2, (3.2) 
Bc{t) = V7 [f^ - {<^x), {t)/2] - iX<jy. (3.3) 

Say we wish to stabilize the pure state with Bloch an- 
gle defined in Eq. (2.6), with r = 1 of course. In 
terms of the ground and excited states, this state is 



^o) = cosyle) 



sin y Iff). 



(3.4) 



Now for this state to be stabilized we must have 

A,it)dt + B,{t)dW{t)'^ \eo) oc \eo). (3.5) 

We cannot say the left-hand-side should equal zero be- 
cause a change in the overall phase still leaves the phys- 
ical state unchanged. However, we can work with this 
equation, and simplify it by dropping all terms propor- 
tional to the identity operator in Ac{t) and Bc{t). We 
can also demand that it be satisfied for the deterministic 
and noise terms separately, because dW(t) can take any 
value. This gives the two equations 

i^a~iXay)\9o) (^leo), (3.6) 
[7 (^—a''a + sin^oo') ~ i2aay 

+ X^{-tsmeoay/2-a^a)] \9o) cx |0o), (3.7) 

where we have put (cTx)^ (t) equal to sin^o, its value for 
the state \6o). 

Solving the first equation easily yields the condition 



A 



-^(1 



(1 + cos 6*0). 



(3.8) 



Substituting this into the second equation gives, after 
some trigonometric manipulation, the second condition 



a 



1 



sin 00 cos ^0- 



(3.9) 



= Xayl{t). 



(3.1) 



These set of driving and feedback amplitude can stabilize 
the system in any particular state in the Bloch sphere (ex- 
cept those on the equator; see below) for unit-detection 
efficiency and Markovian feedback. These features are 
illustrated in Fig. 2 as a full circle (minus the points on 
the equator) which is the locus for the equations with 
optimal driving and feedback. This is contrary to the 
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conclusion of HMH |2J], based on a linearized stability 
analysis, that "long term stability of . . . inverted states 
[i.e. states in the upper half plane] cannot be achieved." 

Once again, it is convenient to construct the stochas- 
tic master equation (SME). This can be d erived from the 
SSE in the same way as before [Eq. The result is 

RP5| 



dpc = dt^'D[a]pc - idta[ay,pc] 

- idtX[ay,apc + Pccr^] + dt^X^ /^)'D[ay]pc 



-dW{t)n[^cr-iXay]pc 



(3.10) 



Also as before, this is an Ito stochastic equation, which 
means that the ensemble average can be found simply by 
dropping the stochastic terms. This time, the result is 
not the original master equation, but rather the feedback- 
modified master equation 



P 



-ilaay, p] +T>[^(T — iX(Ty]p = Cp. (3.11) 



Here we have put the Liouvillian superoperator £ in a 
manifestly Lindblad form. 

It can be shown that the pure state p = \9o){6o\ is a 
stationary solution of this master equation, as expected. 
However, if and only if 9q — ±7r/2, it is not a stable solu- 
tion. That is because the master equation for = ±7r/2 
(it is the same for both cases) has more than one null 
eigenvalue, and any mixture of the two equatorial states 
will be a stationary solution. This means that states on 
the equator, which are equal superpositions of excited 
and ground states, cannot be well-protected against de- 
coherence. We will return to this in Sec. IV C. 



IV. NON-MARKOVIAN FEEDBACK 

A. Non-Markovian feedback SME 

Now we consider the case of non-Markovian feedback. 
This would occur experimentally unless the response of 
the feedback apparatus is flat in frequency space over a 
bandwidth much larger than any relevant system rate. 
For simplicity, we consider the case of just having a time 
delay r. That is, we take the feedback Hamiltonian to 
be 



Hfb = Xayl{t - t). 



(4.1) 



The Stochastic Master Equation (SME ) of the total con- 
ditioned evolution of the system for r finite is given by 

dpc = dt'~iT)[a]pc - idta[(Ty,pc\ + dW{t)y/jH[cr]pc 

- Xdtl{t - T)i[c7y, p^] - dtiX^ /2)[cry,[ay, p^]]. (4.2) 



Here /(s) is still given by Eq. ( |2.14| ), so that the noise 
dW{s) appears twice in Eq. (4.2), with two different time 
arguments, s — t and s = t — t. 



In Eq. ( [l.2| ) it is no longer possible simply to set the 
noise dW equal to its expectation value of zero to obtain a 
deterministic equation. That is because the dW{t—T) ap- 
pearing in the feedback through I{t—T) is correlated with 
the system state, since it already conditioned it through 
the measurement term at time t — t. Also, the atomic 
dynamics are nonlinear, so there is no option but to solve 
Eq. (4.2) numerically. The most convenient way to treat 
the stochastic dynamics is, in general, to use stochastic 
Bloch Equations (SBE). These are simply the stochastic 
equations for the conditioned Bloch vector, defined by 



(/ -f XcCTj, + yc<7y + Zc<7z) . 
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The SBE corresponding to Eq. (|4.2|) are 



(4.3) 





Uc \ dt — dt 



2a + 2XI^r 
-7/2 
-7 - 2A2 






1 



(4.4) 



Here we have used I-r as short-hand for I{t — r). 

We can use these equations to perform the numerical 
simulation of the non-Markovian process. In the numer- 
ical simulation, we have not optimized the driving and 
feedback amplitude to achieve the best possible result for 
a given time delay r. This is as opposed to the Marko- 
vian case with inefficient detection where we optimized 
for every "q. The difference is because in the present case 
the optimization problem is considerable more difficult. 
However, since we expect that the feedback mechanism 
will only work effectively with a small feedback delay 
time, it is thus convenient and reasonable to use the same 
set of driving and feedback amplitude as in the Marko- 
vian case, namely 



A 



-^(1 



-f cos 0o), 



7 

a — — sin ^0 cos ^o- 



B. The effect on stability and purity 



(4.5) 
(4.6) 



Note from Eq. (4.4) that j/c = is a stationary solu- 
tion. Assuming this value for yc, we can again use the 
polar coordinates and r as in Eq. ( ^.6| ). Their equations 
of motion are found from 



+ de = tan" 



dXr 



dzr 



r + dr = \J {xc + dxc)^ + (zc + dzc^ 



(4.7) 
(4.8) 
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Using the Taylor expansion and Ito rules, and assuming 
that the state is initially pure (r = 1) we find 



and 



de 



where 



dr = 



2a + 2^\I[t - r) + ^(2 + cos6')sin6' 



(4.9) 



^{1 + cos e)dW{t), 



dt 



(4.10) 



I{t - T)dt = y/^sin0{t - T)dt + dW{t - r). (4.11) 

That is, the state remains pure and the simulation re- 
duces to t he si ngle non-Markovian stochastic differential 
equation ( I.IC ) which is readily solved in Matlab. 

The two plots in Fig. 3 are the typical quantum trajec- 
tories in Bloch space for 6*0 = 7r/6, starting at the ground 
state. The main plot with more dramatic fluctuations is 
the trajectory via the path through the ground state with 
feedback time delay r = 0.027~^, while the plot with an 
inset into this figure is the one for feedback with no time 
delay and the trajectory is via the path through the ex- 
cited state. Clearly from these two trajectories, there are 
two ways for the system to reach a desired state oiO = 6q. 
As shown in this figure, the trajectory of the main plot 
with r = 0.027~^ continues to evolve stochastically after 
the transients have died away. The parameter 9 wanders 
aromid the desired the pure state = 7r/6 even when 
the system is in steady state. By contrast, that of the in- 
set for Markovian feedback stops precisely at the desired 
pure steady state with no fluctuations persistent. More- 
over, we find that the longer the feedback delay time is, 
the more dramatic the fluctuations are. Therefore, the 
longer time delay will result in more randomness in the 
steady state quantum trajectory. This indicates that the 
stability of the single trajectory is reduced by the effect 
of non-zero feedback delay time. 



Another typical trajectory, this time for 0q = 0, is 
shown in Fig. 4 in Bloch space for t S [0, 1007~^] with a 
feedback delay time r = 0.027~^, again starting at the 
ground state. Firstly, to consider the stability of the sin- 
gle trajectory, we see that the initial evolution is erratic, 
then on a time scale of a few 7^^ the system relaxes 
towards the desired state 6*0 = 0. As in the previous 
case, the fluctuations are still persistent even in steady 
state and there is considerable variation around the de- 
sired angle of 9q = 0. Secondly, to consider the purity 
of the single trajectory, we can see from this plot that 
in a single trajectory the system state is always pure; it 
lives on the surface of the Bloch sphere. However, an 
ensemble (or time) average of the system state will put 
it inside the Bloch sphere. The time-averaged state from 
this trajectory (ignoring transients) is represented by the 
single star dot underneath the point — Q. The larger 
the fluctuations around the closer this point would 
be towards the center of the sphere and the lower the 
purity of the average state. Therefore, we can see that 
the purity of a single trajectory is not be affected by the 
non-Markovian process while the purity of the average 
state is. Rather, the purity of the average state reflects 
the stability of the individual trajectory. 



0.8- 
0.6 
0.4 
0.2 

0- 
-0.2- 
-0.4- 
-0.6- 
-0.8- 



6 1.6- 



8 10 



FIG. 4. Typical single quantum trajectories in Bloch space 
for t G [0, 1007"^] with r = 0.027~^ and 6*0 = 0, starting at 
the ground state. Note that the single star dot underneath 
the single trajectory path indicates the loss of purity of the 
average state, while the single trajectory evolves stochasti- 
cally along the outside full circle which is the locus of pure 
states. 



' ' ■ ' ■ ■ ■' ■ ■ .1 ■ ■.■ 

2 4 6 8 10 

time 

FIG. 3. Typical quantum trajectories for 6 (representing 
the conditioned state), for t £ [0, 107~^] and 9q = 7r/6, start- 
ing at the ground state. Note that the single trajectory can 
go either via the ground state or the excited state to get to 
6(1 — 7r/6. The main plot with a more dramatic fluctuations 
is for feedback with a large time t = 0.027~^, while the inset 
is for feedback with no time delay. 



Fig. 5 is the plot of the locus of the time-averaged 
states in Bloch space for three different feedback delay 
times r = 0, 0.027"-'^, and 0.27~^. For each time delay we 
plot many points, each corresponding to a different value 
of 6q, the angle at which the Markovian feedback would 
stabilize the state. We have deliberately symmetrized 
this plot by combining the simulations of positive and 
negative 6. 



7 



A number of points are worth noting. First, and most 
obviously, the degree of purity (measured by the r, the 
distance from the origin) decreases with r. Second, the 
gap at the equator for r = quickly widens for larger de- 
lay time r, so that the purity of the time-averaged states 
with 6 close to 7r/2 is decreased, which indicates that the 
states near the equator can not be well-protected against 
decoherence. Third, the purity of the time-averaged 
states in the upper half of the Bloch sphere is affected 
much more by the increase of delay time r than those 
in the lower half. Since for very long time delay the 
feedback could not be effective, and since the stationary 
states with no feedback are confined to the lower half of 
the Bloch sphere, this is perhaps not surprising. Fourth, 
the bottom half curves for the larger delay time are closer 
to the ground state even than the no-feedback curve of 
Fig. 2. This is due to the fact that we have not optimized 
the driving and feedback amplitude to search out the best 
possible pure state. With a proper choice of driving and 
feedback amplitude, the curves should always lie between 
the no-feedback ellipse of Fig. 2 and the r = 1 pure state 
circle. 



C. Exception: the equatorial states 



0.5 




As mentioned in Sec. Ill, the story for trying to sta- 
bilize the equatorial states — ±7r/2 is special. First, 
consider the Markovian case. The SBEs for the equato- 
rial states are 



/ dxc \ 


r 





dye 


= 


-7/2 










VldW{t) 





/ Xc 




Vc 


-7/2/ 




l-xl\ 




-XcUc 









(4.12) 



Both Zc and j/c will decay to zero (as required for \6\ = 
7r/2), and their noise terms vanish at that point. By con- 
trast, the equation for Xc is independent of the others, 
and is purely stochastic: 



dx,^^dW{t){l-xl). 



(4.13) 



N 



-0.5 



-1 -0.5 0.5 1 

X 

FIG. 5. Locus of the solutions to the Bloch equations for 
feedback with different values of delay time r. From the out- 
side in, we have r — 0,0.027~^, and 0.2"f~^. The plot with 
r = is an analytical solution, and the others with error bars 
are from numerical simulations. The simulation time for each 
time-averaged state is 10*7"^, giving an effective sample size 
of order 10*, since 7~^ is of order the correlation time of the 
system. 

All of these features (except the last, which was 
avoided through proper optimization) also arose from 
feedback with non-unit-efhciency detection, as discussed 
in our previous paper | |2^ . Both non-Markovian feedback 
and non-unit-efficiency detection decrease the degree of 
stability of the system and thus to limit the capability 
of feedback decoherence control. The difference is that 
with non-unit efficiency detection not only the average, 
but also the individual conditioned state of the system is 
impure. 



Clearly the equatorial pure states with Xc = ±1 are 
stationary solutions to this problem. Also, the system 
will tend to one of these states. But it is also clear 
that Xc has no preference to go to either of these states. 
Hence they are not stable. The ensemble average x is 
unchanging under this evolution. Thus a perturbation 
which moves the state from ccc = 1 to ccc = 1 — e say, 
will result in a proportion e/2 of the states ending up at 
Xc — —1, and a proportion 1 — e/2 ending up at Xc = 1. 

The above discussion all refers to the Markovian case. 
Having a time delay will introduce extra noise in the sys- 
tem, as seen in the preceding section. This is a sort of 
perturbation that will disturb the state of the system. 
Thus we expect that with a time delay, even if the sys- 
tem is close to one of the fixed states Xc — ±1, it will 
not stay there. Rather, we expect that it will at some 
point randomly switch to being close to the other fixed 
point. This is indeed what we observe, as we illustrate 
by showing a typical trajectory for = f in Fig- 6, 
again starting from ground state. On a time scale of a 
few 7~^ the system reaches the region of the no-delay 
fixed point = 1, on one side of the equator. Then at a 
later time, about 57""'^ here, it stochastically hops to the 
other side of the Bloch sphere (going through the ground 
state again), ending up near the other fixed point with 
Xc = -1. 
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FIG. 6. Typical quantum trajectories for 6*0 = tv/2 with 
driving amplitude a = 0, feedback amplitude A = — ^77/2 
and delay time r — 0.027^^, shown by 6 (dots) as functions 
of time t £ [0, 207~^]. This plot shows that the system has no 
preference to go to either equatorial point Xc = 1 {6 — n/2) 
or Xc = -1 {9 = 37r/2). 



D. Analytical solution 

As we mentioned in the introduction, the excited state 
6'o = is the only state where we can get both analytical 
and numerical solutions to the non-Markovian feedback 
process. In this subsection, we will derive the analyti- 
cal solution of the excited state and compare it with the 
numerical one. 

For short time delay we expect that (except for states 
near the equator), 9 will stay near the desired 6*0- Thus, 
it makes sense to define 69 = — 9 o an d do a pe rtur- 
bati ve exp ansion for small 69 in Eq. (|4.5| ), Eq. (4^) and 
Eq. (4.10). This gives the Ito equation 



e exp ansion for small 69 in Eq. (|4.5| 
4.10). This gives the Ito equation 

t) = -7(1 + cos 6*0) cos 9o69{t - t) 
1 



69lt) 



Now the purity of the averagi 
P = 2Tr[p2] - 1 = r 



+7(cos6io + ^cos29o)69{t) ~ y/^sm9o69£_{t) 
-^{l + cos9o)m-T)-at)]- (4.14) 
;e state can be defined as 
' + {ycf + {zcf- (4.15) 



= = (Xr 



Using Uc = 0, Xc = sin 6*, 
(69) = yields, to leading order, 

P^l-{i69f) 



cos 9 and assuming 



(4.16) 



Thus if we can find an expression for {{69)^') we can find 
the purity of the average state. The decrease in purity 
due to fluctuations in 69 can be seen directly in Fig. 4, 
for the case 6*0 = 0. 



Inspecting Eq. ( 4.14 ) we can see that it is solvable only 
if sin 6*0 = so that the multiplicative term 69£^{t) van- 
ishes. This will be the case only when one is trying to 
put the atom in the excited state 9o = 0, or the ground 
state 6*0 = tt. The latter is trivial, as it is achieved exactly 
with no feedback or driving. The former is the interesting 
case, and with = 0, Eq. (4.14) becomes 



r) + ^69{t) 



(4.17) 



We show in the appendix that we can solve this equation 
exactly in Fourier space, and use the result to evaluate 
((50)^) approximately, in the limit r <C 7~^. The result 
is 



P=l- AjT. 



(4.18) 




time delay 

FIG. 7. Plot of the purity against time delay for the atom 
near the excited state. The solid line is the theoretical result 
while the dotted line with error bars is the numerical one. 
The analytical result agrees very well with the numerical one 
when feedback delay time is less than 0.177~^. 

We plot in Fig. 7 the purity against delay time for the 
atom with 9o = 0, according to the analytical approxima- 
tion (4.18) and from numerical simulations. Clearly the 
analytical result agrees very well with the numerical one 
when the feedback delay time is less than about 0.177~^. 
When the feedback delay time is larger than 0.177~^, a 
significant difference appears, and increases as the feed- 
back delay time becomes larger. However, as we can see 
from Eq. (4.18), the feedback delay time must be less 
than 0.257^^ in Eq. ( 1.18 ) in order to satisfy P > 0. 
Therefore for the atom near the excited state, our ana- 
lytical solution is a good approximation for a surprisingly 
large delay time r. 



V. DISCUSSION 

In this paper, we have given a rigorous analysis of an 
anti-decoherence feedback scheme in a two-level atom in 
the limit where the feedback time delay is significant. We 
use numerical simulations to determine the infiuence of 
the non-Markovicity on the effectiveness of the feedback. 
We find that, unlike the case for Markovian feedback, it 
is not possible to stabilize the atom in a fixed pure state 
(except the ground state which is trivially pure by setting 
the driving and feedback to zero). Although the condi- 
tioned system state always remains pure (since we are 
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assuming perfect detection) it does not remain stable but 
continues to wander around the desired fixed pure state. 
Thus the ensemble-averaged state (which is the time- 
averaged state in steady state) is not pure, but mixed. 
As expected, the longer the feedback delay time is, the 
more dramatic the fluctuations are and consequently the 
less pure the average state is. 

We find that the purity of states in the upper half of 
the Bloch sphere is affected much more by feedback delay 
than those in the lower half. This is not unexpected, as 
with no feedback states in upper half cannot be produced 
at all. Also, the purity of states near the equator of the 
Bloch sphere is very much affected. This is because the 
feedback algorithm for stabilizing states on the equator 
cannot distinguish between diametrically opposite points 
on the equator. Because delayed feedback cannot per- 
fectly cancel the measurement back-action noise, trying 
to stabilize at a point on the equator will give a condi- 
tioned state which flips at random from one side to the 
other, on a time scale that decreases as the time delay 
increases. The unconditioned (or time-averaged) state is 
subsequently almost completely mixed. All of these fea- 
tures are qualitatively similar to the results obtained for 
no time delay but with inefficient detection in Ref. [ p2| . 
The major consequence is that states which are near- 
equal superpositions of excited and ground states cannot 
be well-protected against decoherence. 

There is one case in which it is possible to obtain an- 
alytical results, namely when one tries to stabilize the 
atom in the excited state. Our approximate analytical 
result that the purity (the square of the Bloch radius) 
of the average state is given by P = 1 — A'jt (where r 
is the delay and 7 the spontaneous emission rate) agrees 
very well with our numerical results even for quite large 
T, up to 0.177^^. This shows the usefulness of the quan- 
tum trajectory approach to feedback. Not only is it the 
only practical way to treat a non-Markovian problem, it 
can also yield analytical solutions in some cases. Other 
approaches to non-Markovian feedback can only be used 
for systems with linear dynamics, and it is not clear that 
a linearization procedure would work even for the special 
case of stabilizing near the excited state. 

The control of decoherence has become an important 
topic in recent years in many related research areas. In 
quantum computation decoherence is a main limiting fac- 
tor. In a quantum computer information is not stored as 
bits, but rather as qubits [|6|. A qubit is a quantum sys- 
tem with a two dimensional Hilbert space, such as a two- 
level atom. The logical basis states are two orthogonal 
basis states, |0), |1), but the qubit may be in a superposi- 
tion of these states. A quantum computer would consist 
of a large number N of qubits, and could exist as an ar- 
bitrary state in the 2'^-dimensional Hilbert space. If any 
individual qubit undergoes decoherence, such as sponta- 
neous emission, the quantum computation is destroyed. 
Furthermore this decoherence rate scales linearly with 
the number of two-level systems. If quantum computa- 
tion is to become a reality a means must be found to 



control such decoherence. Sophisticated quantum error 
correction methods have been developed to detect and 
control arbitrary amplitude and phase errors |37|. While 
our decoherence-control scheme is by no means as gen- 
eral as error correction, it may be of use for maintaining 
qubits in a particular pure state until required. 

We close by saying a few words about experimental 
realizability. The results of this work shows that a time 
delay r is not fatal to controlling decoherence even in a 
nonlinear system like the two-level atom. The purity of 
the state will, in most cases, remain close to unity as long 
as T is much less than the atomic lifetime 'y~^. This is 
feasible with very fast electronics. The greater difficulty 
is with obtaining high efficiency in detection. Collecting 
a large proportion of the light emitted by an atom in free 
space is very difficult. It is much easier to collect the 
light emitted from a cavity, as this propagates in one di- 
rection. Therefore, the most likely scenario for realizing 
the scheme would be in the context of cavity QED |Q . If 
a two-level atom is strongly coupled (g) to a single cav- 
ity mode, which is strongly damped (k), the combined 
system acts like an effective two-level atom. The output 
beam of the cavity is effectively the spontaneous emis- 
sion of the atom. Then, given that the time delay in the 
feedback loop is relatively small compared to the effective 
atomic lifetime ^ n/g'^), we can control the decoherence 
of the cavity QED system. 
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APPENDIX 

For the atom in the excited state, the SBE in the time 
domain is given by 



-2V7K(i-T)-eW]- 



(6.1) 



Applying the Fourier transformation, we can solve the 
SBE in the frequency domain: 



V7(e'"" - l)e(^) 
[7(1 - 2e*"^) + ico] 



(6.2) 



Applying the inverse Fourier transformation to go back 
to the time domain, we get 



1 f°° 
d0{t) = — / dLue 



-..t V2(f!l^L^(^). (6.3) 
7(|-2e''^^)-Fia;^^ ' ^ ' 
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Since (^(lu)(^{uj )) — 2TrS{Lu + uj ), we find 



{SO' it)) = 



47(2 — 2 cosujT)duj 



7^(1 — 2 cosuJTy + {uj — 27sinwT)2 



16 sin^ {ujt/ 2) duj 



1 

2^ 
1 

'^1 Jo (| — 2 coswt)^ + (07/7 — 2 sinLJT)2 

The above integral is too complicated to be useful. 
Also, we are really only interested in the limit r ^ 7""'^, 
since that is the limit in which SO is expected to be small 
so that the approximations leading to Eq. (^^) will be 
valid. With this in mind, consider a parameter c such 
that 7 ^ c ^ r^^, and c^r <C 7^. Then spht the uj- 
integration up as = + . When a; < c, we have 
16sin^ ujt/2 < (2a'T)^<Cl in the numerator, while the de- 
nominator is greater than or equal to j. Therefore the 
integration from to c is bounded above by IGt^c^ /3. On 
the other hand, when w > c, we have w ^ 7 so that w/7 
is much greater than all the other terms in the denomi- 
nator. Therefore we can approximate the integration in 
this region as 



duj 



16Y sin^(wT/2) 



(6.4) 



which, since c <C t ^, can be further approximated as 
1672sin2 (wt/2) 



duj- 



(6.5) 



Since by assumption c'^r <C 7^, this integral is always 
much greater than that from to c. We thus finally ar- 
rive at 

((56'^(t)) ~ — X 7247rT = 47r < 1. (6.6) 
7r7 
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